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PREFACE 

This  Memorandum  results  from  RAiro' s  continuing 
study  of  computational  aspects  of  the  solution  of  the 
equations  of  mathematical  physics.  A  new  basic  partial 
differential  equation  of  particle  mechanics  is  derived. 
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SUMARY 

Consider  the  motion  of  a  particle  on  a  line  where ^ 
as  usualj  vie  characterize  the  process  by  a  Hamiltonian 
H  =  H(q,p).  The  equations  of  motion  are 

(1)  q  =  Hp. 

(2)  -  P  =  0  <  t  <  T. 

As  boundary  conditions  we  specify 

(3)  <l(0)  =  0, 

(4)  p{T)  =  c. 

If  our  aim  is  to  solve  the  system  of  equations  (l)  and  (2), 
subject  to  the  boundary  conditions  in  (3)  and  (4),  using  a 
digital  computer,  we  face  the  well-known  difficulties  of 
solving  nonlinear  two-point  boundary— value  problems. 

It  is  tempting  to  try  to  determine  the  unlcnown  dis¬ 
placement  at  time  T,  q(T),  so  that  we  shall  have  a 
complete  set  of  conditions  at  time  T.  This  would  enable 
us  to  integrate  the  system  (l)  and  (2)  numerically, 
subject  to  initial  values.  Toward  this  end  we  introduce 
the  function 

(5)  r(c,T}  =  the  displacement  at  time  T,  the 

initial  displacement  being  zero  and 
the  momentum  at  time  T  being  c. 


We  derive  a  partial  differential  equation  for  this  function. 


f 


-vil- 

CONTEKT?S 


PREFACE . Ill 

SUM-IARY . V 

Section 

1.  INTRODUCTION  .  1 

2.  INVARIANT  IMBEDDING . 2 

3.  AN  EXAMPLE . 3 

4.  DISCUSSION . 4 

REFEPcENCES . 5 


A  NOTE  ON  HAIvlILTON'S  EQUATIONS  Al^D  INVARIANT  IlffiEDDING 


1 .  INTRODUCTION 

Consider  the  motion  of  a  particle  on  a  line  where, 
as  usual,  we  characterize  the  process  by  a  Hamiltonian 
H  =  H(q,p).  The  equations  of  motion  are 

(1.1)  q  =  Hp, 

(1.2)  -  P  =  H  -  0  <  t  <  T. 

As  boundary  conditions  v;e  specify 

(1.3)  q(0)  =  0, 

(1.4)  p(T)  =  c. 

If  our  a.im  is  to  solve  the  system  of  equations  (l.l)  and 
(1.2),  subject  to  the  boundary  conditions  in  (1.3)  and 

(1.4) ,  using  a  digital  computer,  we  face  the  well-kno;'/n 
difficulties  of  solving  nonlinear  two-point  boundary- 
value  problems. 

It  is  tempting  to  try  to  determine  the  unknown 
displacement  at  time  T,  q(T),  so  that  we  shall  have  a 
complete  set  of  conditions  at  time  T.  This  would 
enable  us  to  integrate  the  system  (l.l)  and  (1.2) 
numerically,  subject  to  initial  values.  Toward  this 
end  we  introduce  the  function 
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(1.5)  r(c,T)  =  the  displacement  at  time  T,  the 

initial  displacement  being  zero  and 
the  momentmii  at  time  T  being  c. 

We  wish  to  derive  an  equation  for  r(c,T). 

2 ■  INVARIANT  IMBEDDING 

In  earlier  papers  [l^2],  we  showed  the  following: 
If 


(2.1) 

^  =  F(u,v),  u(0)  =  0 

(2.2) 

-  =  0(u,v),  v(x)  =  0 

for 

(2.3) 

0  <  z  <  X, 

and  if 

(2A) 

r(c,x)  =  u(z)|  , 

z=x 


then  the  function  r(c,x)  satisfies  the  partial 
differential  equation 

(2.5)  r^  =  F(r,c)  +  G(r,c)r^ 

and  the  initial  condition 


(2.6) 


r(c,0)  =  0. 
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Upon  applying  this  result  to  the  system  of  equations 
in  (l.l)  and  (1.2),  we  find  that  the  function  r(c,T) 
satisfies  the  equation 

(2.7)  ^  =  Hp(r,c)  + 

or,  more  elegantly, 

(2.8)  = 

This  is  the  desired  result,  apparently  a  new  equation  of 
mechanics. 

3.  AN  EXAMPLE 

Consider  the  case  of  harmonic  oscillations  for 

which 

2 

(3*1)  H(q,p)  =  ^  +  ^  kq^. 


and 

(3.2)  q(0)  =  0,  r(T)  =  c. 


For  the  unknown  displacement  at  time  T,  r(c,T)j 
equation  (2.8)  becomes 


(3.3) 


Sr  _  S  fc^  1 
BT  -  L2i  +  2 


=  5  +  *''’ 


Sr 

Be’’ 


The  solution  of  this  equation,  subject  to  the  condition 
(2.6)  is 


This  provides  the  desired  displacement  at  time  T  and 
also  shows  that  the  function  r(c,T)  may  become  Infinite 
for  a  finite  value  of  T. 

4 .  DISCUSSION 

The  previous  discussion  may  be  generalized  to  the 
case  of  many  particles  moving  in  three-dimensional  space 
using  the  theorem  in  111.  Furthermore,  various 
combinations  of  the  q's  and  p's  may  be  specified  at  the 
ends.  There  are  immediate  applications  to  perturbation 
analysis  of  results  of  this  nature;  see  [  3^  . 
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